Superlattice Patterns in the Complex Ginzburg-Landau 
Equation with Multi-Resonant Forcing 

Jessica M. Conway 1 and Hermann Riecke 1 ' 2 
March 3, 2008 

1 Engineering Sciences and Applied Mathematics, Northwestern University, Evanston, IL 60208, 
USA 

2 Northwestern Institute on Complex Systems, Northwestern University, Evanston, IL 60208, USA 

Abstract 

Motivated by the rich variety of complex patterns observed on the surface of fluid layers that 
are vibrated at multiple frequencies, we investigate the effect of such resonant forcing on systems 
undergoing a Hopf bifurcation to spatially homogeneous oscillations. We use an extension of 
the complex Ginzburg-Landau equation that systematically captures weak forcing functions 
with a spectrum consisting of frequencies close to the 1:1-, the 1:2-, and the l:3-resonance. 
By slowly modulating the amplitude of the 1:2- forcing component we render the bifurcation 
to subharmonic patterns supercritical despite the quadratic interaction introduced by the 1:3- 
forcing. Our weakly nonlinear analysis shows that quite generally the forcing function can be 
tuned such that resonant triad interactions with weakly damped modes stabilize subharmonic 
patterns comprised of four or five Fourier modes, which are similar to quasi-patterns with 4-fold 
and 5-fold rotational symmetry, respectively. Using direct simulations of the extended complex 
Ginzburg-Landau equation we confirm our weakly nonlinear analysis. In simulations domains 
of different complex patterns compete with each other on a slow time scale. As expected from 
energy arguments, with increasing strength of the triad interaction the more complex patterns 
eventually win out against the simpler patterns. We characterize these ordering dynamics using 
the spectral entropy of the patterns. For system parameters reported for experiments on the 
oscillatory Belousov-Zhabotinsky reaction we explicitly show that the forcing parameters can 
be tuned such that 4- mode patterns are the preferred patterns. 

1 Introduction 

Complex but ordered spatio-temporal patterns, characterized by multiple length scales, have been 
observed and studied in a range of systems. In particular in the Faraday system, in which a thin 
layer of fluid is vertically vibrated leading to patterns on the surface of the fluid, various kinds of 
superlattice patterns and quasi-patterns have been found experimentally depending on the frequency 
content of the forcing function [1-9]. Related complex patterns have also been observed in optical 
systems [10], in vertically vibrated fluid convection, where they arise from the competition of different 
instability mechanisms [11], and on the surface of ferrofluids driven by time-periodic magnetic fields, 
where they are due to spatial period-doubling [12]. 

From a purely spatial point of view the quasi-patterns as well as the superlattice patterns observed 
in the above forced dissipative systems are closely related to quasi-crystals obtained in thermody- 
namic equilibrium [13]. There the significance of damped, resonating modes for the stabilization 
of quasi-crystals has been recognized [14,15]. Recently, two-dimensional quasi-crystals have been 
found in soft-matter systems [16]; their dodecagonal symmetry has been discussed in the context of 
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a Swift-Hohenberg-type model [17] that had previously been introduced to get insight into Faraday 
patterns [18]. 

The temporal aspect introduced by the forcing of non-equilibrium systems greatly increases the 
richness of patterns observed. To wit, contributing to the great variety of complex patterns observed 
in the Faraday system is the ability to control in detail the temporal wave form of the vibration's 
forcing function [3,5,7,19,20]. This allows extensive tuning of the interaction between plane waves of 
different orientation, which can in turn stabilize superlattice patterns and quasi-patterns. Two main 
stabilization mechanisms have been identified, both of which exploit the presence of weakly damped 
resonating modes. In one case the competition between plane-wave modes of different orientation is 
suppressed for a relatively narrow range in the angle subtended by the competing modes [19-22,53]. 
In the other case the self-coupling of each mode is strongly enhanced, rendering the competitive 
coupling between the modes effectively weak over a considerable range in the angle [22-28] and 
allowing patterns comprised of multiple modes to become stable [29] . The latter mechanism captures 
qualitatively the scenario envisioned for the stabilization of 'turbulent crystals' [15]. We will focus 
on this second mechanism and will make use of the fact that the temporal forcing introduces an 
additional symmetry that can be exploited to select or suppress certain types of patterns. 

In this paper we investigate, motivated by the richness of patterns observed in the Faraday 
system, whether such complex spatio-temporal patterns are also accessible in systems undergoing a 
Hopf bifurcation to spatially homogeneous oscillations. These systems constitute a class that differs 
from those in which quasi-patterns have been observed previously and they may offer the potential 
for additional complexity through the interaction between the spontaneous oscillation arising from 
the Hopf bifurcation and the external forcing. Chemical oscillations like those observed in the 
Belousov-Zhabotinsky reaction are a classic example of such a system. 

In the absence of any temporal forcing the Hopf bifurcations we have in mind lead to spatially 
homogeneous oscillations or long- wave traveling waves, which may break up to form spirals or more 
complex chaotic states (e.g. [30-34]). With temporal forcing, which in the case of the Belousov- 
Zhabotinsky reaction can be achieved by time-dependent illumination to exploit the photosensitivity 
of the reaction, the oscillations can become locked to the forcing and - depending on the ratio 
between the Hopf frequency and the forcing frequency - different types of patterns can arise [35-40] . 
Thus, forcing near twice the Hopf frequency can lead to competition between two types of domains 
differing in their temporal phase or to labyrinthine patterns. For a frequency ratio of 1:3, labyrinthine 
patterns, spirals, or competing domains have been observed. Ordered patterns with multiple length 
scales have been observed experimentally in chemical systems so far only when a spatially periodic 
illumination mask was applied to initialize the pattern [41,42]. We focus on the case of spatially 
uniform illumination. 

The organization of the paper is as follows. In Sec(2]we extend the complex Ginzburg-Landau 
equation (CGLE) to include the terms that describe the external forcing at various frequencies. The 
derivation of this equation for the Brusselator, which is a simple model for chemical oscillations, is 
sketched in the Appendix. In Sec|3]we present a linear stability analysis of the equation to find the 
onset of standing waves that are phase-locked to the driving. In particular we focus on the vicinity 
of the codimension-2 point at which the subharmonic and the harmonic standing waves bifurcate 
simultaneously. To determine the stability of the desired subharmonic patterns we derive in Sec(4] 
the corresponding amplitude equations by performing a weakly nonlinear analysis of the CGLE, and 
then use energy arguments to guide us in terms of the relative stability of various pattern comprised 
of different numbers of modes. In order to confirm our predictions for the pattern selection, we 
perform in Sec (5] numerical simulations of the CGLE in small and large domains, and characterize 
the temporal evolution of patterns by using a spectral pattern entropy. A brief account of these 
results has been published previously in [43]. 

2 The Complex Ginzburg-Landau Equation 

We are interested in the formation of complex patterns in systems that undergo a Hopf bifurcation 
at vanishing wave number and that are forced at frequencies near integer multiples of the Hopf 
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frequency. Near the Hopf bifurcation and for weak forcing such systems can be described by a 
suitably extended complex Ginzburg-Landau equation (CGLE), the form of which can be derived 
using symmetry arguments. 

In the absence of forcing the complex amplitude A(t) satisfies the usual CGLE, 

OA 

— =a 2 A + a 3 AA + a 4 A\A\ 2 , (1) 

where a$ G C, i = 1...4 and A is the Laplacian in two dimensions [33]. The CGLE exhibits 
the normal-form symmetry T T : A — > Ae luJT reflecting the invariance of the original system under 
translations of the time t by an arbitrary amount r, T T : t — > t + r. Here uu is the Hopf frequency. 
The slow time t is given by t = J 2 t, 0<(S<1. 

The extension of the CGLE that describes strongly resonant multi- frequency forcing of the form 
F = fie ZUJt + f 2 e 2l,UJt + fse 3zuJt + c.c. can be obtained by considering the analogous center-manifold 
reduction of the extended dynamical system in which the forcing amplitudes /i, f 2l and fs are 
considered as dynamical variables that vary on the slow time scale t. Under time translations T T 
they transform as f\ — > fie %UJT , f 2 — > f2e 2zuJT , fs — > f^ %UJT (e.g. [44]). To cubic order in A the most 
general equation that is equivariant under T r is then given by 

— = ai + a 2 A + a 3 AA + a 4 A|A| 2 + a 5 A + a 6 A 2 , (2) 
at 

where ai = b X \f\ + 612/2/3, «2 = 621 + ^22 1 /3 1 2 , «5 = 651/2, «6 = &6i/3- The 6^- are 0(1) complex 
coefficients. For all terms to appear at the same order in the final amplitude equations, we use the 
scaling A(t) = 0(S), f\ — 0(J 3 ), f 2 = 0(5 2 ), and fs = 0(6). Spatial coordinates are scaled as 
(x,y) = 5(x,y), where (x,y) are the original spatial coordinates of the system. 

The forcing terms fj satisfy decoupled evolution equations on their own. In the simplest case 
this evolution expresses a detuning Vj of the forcing fj from the respective resonance and the fj 
satisfy 



dfj 



ivjfj, j = 1...3. (3) 



In general, the detuning introduces time dependence into 

In a previous study [45] we investigated the simpler case in which the forcing parameter f\ 
satisfies 611/1 = —612/2/3 and the detuning parameter ^3 satisfies ^3 = 3^/2, which allowed us to 
make all coefficients of the extended CGLE (I2| time-independent by a transformation A — > Ae tU2t / 2 . 
For that case we found that the possibly stable superhexagon pattern, a superlattice pattern, arises 
in a transcritical bifurcation. We showed that the amplitudes on the upper branch of the transcritical 
bifurcation are 0(1) and that therefore a weakly nonlinear analysis taken to cubic order is insufficient 
to describe these patterns. The bifurcation is transcritical because the quadratic term in the CGLE 
introduces also a quadratic term in the amplitude equations describing the competing Fourier modes. 
To eliminate the latter quadratic term without losing the resonant triad interaction we choose in 
the following the forcing near twice the Hopf frequency to be quasi-periodic, 

p p 1 p ., u df 2 i . j. df 22 . » / a \ 

H hi + 722 With ZI/21/21, ^22/22- (4) 

The difference between the two detunings v 2 \ and v 22 introduces then a periodic time dependence 
in ([2]). To obtain the desired patterns, we exploit the spatiotemporal resonances induced by the 
time dependence and focus on patterns that are subharmonic in time. The quadratic interaction of 
two subharmonic modes induces a harmonic mode and does not generate a quadratic term in the 
amplitude equation for the subharmonic mode. Consequently, we expect a pitch-fork bifurcation. 

Considering the quasi-periodic forcing Q we simplify p| by absorbing the time dependence of 
/21, e lU2lt , into A through the transformation A — > Ae lU2 ^ 2 . Further, we write as = 1 + ifl by 
rescaling the spatial coordinates and — — (1 + ia) by rescaling A. We focus on the case of a 
supercritical Hopf bifurcation and choose the real part of to be negative. We now introduce 
restrictions on the forcing fj for the purpose of making the analysis manageable at the expense of 
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some generality. To eliminate the homogeneous term a\ we choose bufi = — &i2(/2i + /22)/3- To 
remove the time dependence from the A 2 term we choose v% = 3^21 /2. This yields 

B A 

— = (,1 + ia)A + (1 + i/3)AA - (1 + ia)A|A| 2 + ( 7l + l2 e l ^- v ^ l )A + 77A 2 . (5) 

Here the growth rate /i is modified by the forcing function I/3I 2 , cr is a linear function of the detuning 
V21, and r] is a complex linear function of We rewrite 77 in magnitude and phase, 77 = pe^, and 
write 71 = 7cos(x) and 72 = 7sin(x) with x characterizing the relative forcing strengths. Finally 
we set v = ^22 — ^21, to get the version of the complex Ginzburg-Landau equation (CGLE) that we 
will investigate in the following, 

3 A 

— = (fjL + ia) A + (1 + i(3)AA - (1 + ia)A\A\ 2 + 7 (cos(x) + sin(x)e^) A + pe^A 2 . (6) 

In this paper we stay below the Hopf bifurcation and focus on the case \i < 0. As in the Faraday 
system, in the absence of forcing, the only solution is then the basic state A = 0. 



3 Linear Stability Analysis 

To obtain the onset of standing waves that are phase-locked to the driving we consider the stability 
of the solution A = of ([6| [36]. Linearizing ([6| about A = and splitting the equations into real 
and imaginary parts by setting A = A r + iAi, we obtain linear partial differential equations for A r 
and Ai. Since these equations are time-periodic with period 2tt/u we use Floquet theory in the same 
way as described in [46]. The solutions to the linear problem in Floquet form are 

A \ ~ °° / Y 



/ Ti— — OO x 

Here e^ + ^ 4 is the Floquet multiplier. We restrict £ to two values corresponding to a harmonic 
response (£ = 0) or to a subharmonic response (£ = 1/2). To obtain the neutral curve for the 
standing waves, the ansatz ^ is substituted into the linearized system with (5 = 0, yielding the 
infinite-dimensional system of equations 

/c 2 +7(/c)cos(x)-i(n + (> -a + (3k 2 \ ( X n 



——00 ^ ^ 

, l(k) sin(x) 



a — (3k 2 11 - k 2 — j(k) cos(x) — i(n + Cj v J \ Y n 

(8) 

1 * \( x n+1 
i 1 )\ Y n+l 



2 

which can be written as a matrix equation: 

M(C)v = 7 (fc)Lv, (9) 

where 

v = (... X n _i Y n -\ X n Y n X n+ i 

Here M is block-diagonal and L couples adjacent modes (X n ,Y n ) and (X n ±i,Y n ±i). Equation ([9| 
is a generalized eigenvalue problem for the eigenvalue j(k), which represents the forcing amplitude 
at the onset of the instability for a given k. In our formulation of the linear problem 7 is assumed 
to be real; complex eigenvalues therefore do not correspond to solutions of the original problem. To 
solve for the eigenvalues for a given value of the wave number /c, we truncate the sum in (|9| at some 
N and calculate the eigenvalues. We test for convergence by requiring that for the real eigenvalues 
7^ the total change in all eigenvalues when N is increased by 1 is smaller than a tolerance A, 
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Figure 1: Harmonic (green, diamonds) and subharmonic (blue, circles) neutral stability curves for 
/i = — 1, cr = 4, (3 = 3. The parameters \ an d v are tuned to obtain the desired K. (a) For 
K = 1.857, x = 0.4883189 and v = 3.824884. (b) For K = 2, x = 0.4767180 and v = 4.2. 



XX 7i + —7} ) 2 < A = 10 -10 . We find that generally TV = 10 is sufficient. Repeating this process 
over a range of k we construct the neutral stability curves ^ H \k) for the harmonic mode (( = 0) 
and 7^(/c) for the subharmonic mode (( = 1/2). The global minimum of these curves yields the 
respective critical values (ki H \^i H ^) and (ki S \ 7^). 

The weakly nonlinear analysis below shows that weakly damped harmonic modes have a strong 
impact on the selection of subharmonic patterns via resonant triad interactions. Our aim is to 
exploit this sensitivity to stabilize complex patterns like superlattice patterns. We therefore focus 
here on parameters for which the minimum of the harmonic mode is only slightly above that of 

( H) ( S) 

the subharmonic mode; we use 7c — 7c =0.01. The resonant triads rely on spatiotemporal 

( H) ( S) 

resonances that depend decisively on the wavenumber ratio K = k c /k c • We consider the two 
cases K = 2 cos(tan _1 (2/5)) = 1.857 and K = 2 and show in Fig II] the corresponding neutral 
stability curves for 7^ — 7^ = 0.01. 

Varying the parameter x characterizing the ratio of forcing strengths near the l:2-resonance shifts 

( H) ( S) 

the critical forcing parameters 7 C and j c relative to each other, while varying the parameter v 

( H) ( S) 

characterizing the detuning difference shifts the critical wavenumbers kc and k c . Requiring 

( H) ( S) 

7c 7o 0-01 for a fixed value of K defines a codimension-2 point that determines both x an d 

v. For fixed K = 2 their dependence on the dispersion (3 and the damping \i is illustrated in Fig(2] 

(S) (s) - 

along with the corresponding values of 7 C and k c . Since A satisfies the CGLE with the opposite 
signs of the imaginary parts of the equation, only positive values of (3 need to be considered. 

For x — the critical wave number is given by k 2 c = (fi + v(3)/{l + (3 2 ) corresponding to a 
critical forcing 7^ = (y — n(3) 2 /(l + (3 2 ) [36]. Thus, k c has a local maximum and k 2 c ~ 1/(3 for large 
(3] qualitatively the same behaviour is found at the codimension-2 point, as illustrated in Figj2]i. 
Figj2^ confirms that the amount of forcing necessary to generate patterns increases as one goes 
further below the Hopf bifurcation, i.e. as ji becomes more negative. 



4 Weakly Nonlinear Analysis 

To determine the stability of the desired subharmonic patterns we derive the corresponding ampli- 
tude equations by performing a weakly nonlinear analysis of (jHJ. We then use energy arguments to 
guide us in terms of the relative stability of the various pattern comprised of different numbers of 
modes. We are in particular interested in the stability of patterns comprising three or more modes. 



4.1 Amplitude Equations 

Splitting the CGLE (|6| again into real and imaginary parts, we obtain a system of real partial 
differential equations. For a subharmonic pattern with N modes, we expand {A r ,Ai) about (0,0) 
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Figure 2: Linear parameters at the co dimension- 2 point K = 2 and 7, 
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0.01 for detuning 



a = 4 as a function of the dispersion /3 and the damping /i. (a) Detuning (b) forcing ratio x-> ( c ) 
critical forcing strength 7c , and (d) critical wave number kc 
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(10) 



where < e <C 1 and the complex amplitudes Z m (T) depend on the slow time T 



The 



corresponding wave vectors are denoted by k m . We also expand 7 as 7 = 7 C + e 72. 

A standard but lengthy calculation yields amplitude equations for the Zj(T) describing the 
N— mode pattern: 



dZj 
dT 



A72^i 



bo\Zi\ 



N 

■ E 



b(6 k 



Z u i = l...N, 



(11) 



where Oij corresponds to the angle between the wave vectors Iq and k.-. Note that there are no 



quadratic terms in (11), because these amplitude equations correspond to subharmonic patterns 
and must therefore be equivariant under the symmetry Zj — > —Zj. In order to have stable complex 
patterns made up of many modes, the competition between modes needs to be sufficiently weak. 



Since in the amplitude equations (11) only pairwise mode interactions arise, a minimal condition for 



complex patterns to be stable is the stability of rectangle patterns, which are comprised of 2 modes 
at some angle with respect to stripes. This requires \b(0)/bo\ < 1. 
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The self- and cross-coupling coefficients, bo and b(9ij), respectively, can be strongly influenced by 
resonant triad interactions as illustrated in Figj3| There the black circles represent the subharmonic 
critical circle and the red circle the wavevector of the most weakly damped harmonic mode for a 
given choice of critical forcing 7. At quadratic order resonant triad interaction takes place through 
two mechanisms, each feeding into the coefficients at cubic order. The first mechanism is through 
the interaction of two different wave vectors (ki + k2 in Figj3^i) seperated by an angle r : the 
cross-coupling coefficient b{0) is strongly impacted near = r if the mode excited at quadratic 
order through the interaction of these two wave vectors is weakly damped [19-22]. The angle 6 r is 
a function of the wavenumber ratio K, r = 2 cos -1 (if /2). The second mechanism is through the 
interaction of a wave vector with itself (2ki in Figjsj}): in this case, if the harmonic mode at 2k c is 
weakly damped, the self-coupling coefficient bo is strongly influenced [22,23,25]. In Figj3^i K < 2 
(0 r 7^ 0), so the mode excited by the vector ki + k 2 is weakly damped and the mode excited by the 
2ki vector is strongly damped. In Figj3b K = 2 (6 r = 0), so the opposite is true. 




(a) 




(S) 

Figure 3: Resonant triad interaction. Inner circle represents critical circle at kc , outer circle 
represents weakly damped mode at (a) K = 2tan _2 (2/5), (b) K = 2. The wave vector ki + k 2 is 



relevant for b(6) while 2ki affects bo (cf. eq.(ll)). 



Fig (4] shows the effect of the resonant triad interactions on the self-coupling coefficient bo, by 
plotting bo as a function of the 1:3 forcing phase $ for different values of the 1:3 forcing strength 
p. FigjZ^i corresponds to the case illustrated in Figj3^i with K = 2 cos(tan _1 (2/5)); with this choice 
of K the resonance angle r = 2 tan -1 (2/5) is near 7r/4 and the associated Fourier modes fall on 
a regular grid [47] for later comparison with numerical simulation. In Figj4j} K = 2 and the effect 
of the resonant triad on bo is greatly enhanced - note the difference in scale compared to Figj4^i - 
and even quite weak forcing p can change the sign of bo and with it the direction of the bifurcation 
to stripes. In the following we restrict ourselves to values of the phase <1> for which bo > and the 
stripes bifurcate supercritically. 

To assess the mode competition for K = 2, FigjF] shows the dependence of the cross-coupling 
coefficient ratio b(0)/bo on 0. To enhance the stability of patterns comprised of multiple modes 
by maximizing the range of weak competition we maximize the damping bo and choose $ = 37r/4 
(cf. Figj4|. Since K = 2, the self-coupling bo is enhanced by the resonant triad, while away from 
= r = the cross-coupling b(6) is only weakly affected. Thus, as expected, for p > 1 the ratio 
b(0)/bo is strongly reduced away from = allowing for rectangle patterns corresponding to angles 
as small as = C ~ 30° and smaller to be stable with respect to stripes. Consequently, patterns 



comprised of four or possibly even more modes are expected to be stable, as discussed in Sec. 4.2 



The mode interaction also depends on the linear dispersion a and the nonlinear dispersion pa- 
rameter p. These two parameters are system parameters that may not be tunable in experiments. 
It is therefore of interest to assess over what range in these parameters complex patterns are to be 
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(a) 
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0° 60° 120° 180° 

1:3 Forcing Phase 4> 



Figure 4: Self-coupling coefficient bo plotted against the 1:3 forcing phase $ for various values of the 
l:3-forcing strength p with a = — 1. (a) K = 2 cos(tan _1 (2/5), other parameters as in Figflji; (b) 
K = 2, other parameters as in Fig(T|3. 




30° 60° 90 

Angle Between Modes 



Figure 5: Coupling coefficient ratio b(0)/bo for K = 2 with <I> = 37r/4 with linear parameters as 
in FigfTJ) and nonlinear dispersion a = —1 for different strengths of 1:3- forcing p. Rectangles are 
stable for |&(0)/6 O | < 1- 

expected and to what extent the forcing parameters can be adjusted so as to reach the regions of 
interest for given system parameters. 

The effect of the nonlinear dispersion parameter a on the ratio b(0)/bo is shown in Figj6^i for 
p = 1. Somewhat similar to the dependence on p, with increasing a the ratio b{0)/bo decreases 
and the range of for which rectangle patterns are stable increases. For more positive values of a 
even small forcing amplitudes p reduce the mode interaction significantly, but as p is increased the 
ratio b(0)/bo saturates at values that are not much smaller than shown in Figj5]for a = —1. In the 
simulations discussed in Sec(5]below we use a moderate value of a — — 1. For single-frequency forcing 
it was found that the self-coupling coefficient bo changes sign for a = /?, rendering the bifurcation 
subcritical for a > (3 [36]. Interestingly, with multi- frequency forcing the sign of bo depends also on 
the forcing strength p and, in fact, even for a = (3 the bifurcation to stripes can be made supercritical 
by increasing the forcing strength. This is illustrated in Figj6}} where bo is shown as a function of p 
for different values of a near (3 = 3. 

To assess the dependence of the mode interaction on the linear dispersion coefficient (3 we do not 
perform a scan in (3 but rather focus on one other value, (3 = 1.4, which corresponds to the value 
found experimentally for the Belousov-Zhabotinsky reaction [48]. We adjust the forcing parameters 

( H) ( S) 

X and v to stay at the codimension-2 point K = 2 and 7c — 7q = 0.01. The nonlinear dispersion 
coefficient was found in the experiments to be a = —0.4. Since the functional form of the angle 
dependence of b{0) does not change much with the parameters we use the minimal angle C for which 
rectangles are still stable with respect to stripes as a proxy for the strength of the mode competition. 
Thus b(0 c )/b o = 1. In Fig(7| we show 6 C as a function of the forcing strength p for various values of 
the overall detuning a. Clearly, even for these experimentally relevant parameters the resonant triad 
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Figure 6: Coupling coefficient ratio b(0)/bo (a) and self-coupling bo (b) for K = 2 and = 37r/4 
with parameters as in Fig [I}} for 1:3 forcing strength p = 1 and different values of the nonlinear 
dispersion parameter a. 



interaction can be exploited to make the mode competition sufficiently weak to stabilize rectangle 
patterns over a wide range of angles. This suggests that multi-mode patterns should become stable 
for moderate values of the forcing parameters, which is confirmed below. 



•-•g= 
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Stable Rectangles 



1 1.5 2 2.5 
1:3 Forcing Strength p 



Figure 7: Belousov-Zhabotinsky case. Minimal angle C for which rectangles are stable to stripes 
(b(0 c )/bo = 1) as a function of the l:3-forcing strength p for different values of the detuning a. 
Parameters p = — 1, a = —0.4, /3 = 1.4 with x an d ^ tuned to the codimension-2 point K = 2 and 

(#) (5) n m ' 
7c - 7o = 0.01. . 

For K = 2 cos(tan _1 (2/5)) the resonant triad affects mostly b(6) rather than 6 - Figj8]shows the 
resulting ratio b(0)/bo as a function of the angle 0. To ensure bo > we use again = 37r/4. As 
expected b(0)/bo exhibits a prominent peak near = r = 2 tan -1 (2/5). As in Fig{5]the enhancement 
is increased with p (Figj8^i), and as a becomes more positive (Figj^)). We also see from these 
plots that the range for which b(0)/bo < 1, that is, the range of stable rectangles, is limited to 
60° < < 90°. Thus, we anticipate only stripe, square or hexagon patterns to be stable in this 
regime. 

For K = 2 cos(tan _1 (2/5)) the competition between modes subtending an angle near r is 
strongly enhanced for values of the 1:3- forcing phase near $ = 37r/4 resulting in a suppression 
of # r — rectangles. For the same K but $ = 7r/4, we expect correspondingly a selective enhancement 
of the # r — rectangles as long as p is sufficiently small so that bo > (see FigQ. This is shown in 
Figj9j which depicts the behaviour of the cubic coupling coefficient ratio as a function of the angle 
for different values of the l:3-forcing strength p and of the nonlinear dispersion a in this regime. 
As the 1:3- forcing strength p is increased, and as the nonlinear dispersion a becomes more positive, 
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Figure 8: Coupling coefficient ratios for K = 2 cos(tan _1 (2/5)) and = 37r/4 with parameters as 
in Figflji for (a) nonlinear dispersion parameter a = — 1 and different strengths of 1:3- forcing p and 
(b) 1:3 forcing strength p = 1 and different values of the nonlinear dispersion parameter a. 



the dip with minimum at = r — 2 tan -1 (2/5) becomes more and more pronounced. In parallel 
the range in 9 of stable rectangle patterns increases around r . Due to the narrowness of the dip 
we expect at most stripe and # r -rectangle patterns to be stable. For values of p and a such that 
b(0)/bo < — 1 the corresponding rectangles pattern become subcritical and our weakly nonlinear 
analysis to cubic order is insufficient to make predictions about pattern selection involving those 
modes. 




Figure 9: Coupling coefficient ratios for K = 2 cos(tan _1 (2/5)) and $ = 7r/4 with parameters as in 
Figjl^i for (a) nonlinear dispersion parameter a = — 1 and different strengths of 1:3- forcing p and 
(b) 1:3 forcing strength p = 1 and different values of the nonlinear dispersion parameter a. 



4.2 Competition between complex patterns 

If multiple patterns are simultaneously linearly stable they can coexist and compete in sufflcently 
large systems. Typically they form then domains and the competition involves the motion of walls 
separating the domains. If the system allows a Lyapunov functional this competition can be char- 
acterized in terms of the difference between the energies associated with the respective patterns. 
Unless the domain walls become pinned by the underlying pattern [49,50], it is expected that the 
final state arising from random initial conditions consists of the pattern with minimal energy [28,51]. 
Because = b(6ji), (11) can indeed be derived from a Lyapunov functional such that 
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dT 
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for a pattern with N modes. The Lyapunov function gives the energy Ttf for an equal- amplitude 
N— mode pattern, which we rescale to obtain J 7 ^, 



N 



N 



(7 -7c) 2 b + E 



N 

n=l,n^j 



i b(0 jn ) 



(13) 



for some j G [l,iV]. The sum in the denominator represents the sum of all cubic coefficients in (11). 
Since the modes are evenly spaced, the choice of j in (13) is arbitrary. We use the energy Tj^ as a 



guide to predict which N- mode pattern will ultimately emerge. 



Figs fTolli 



show how the energies Tjsf vary with p, the strength of forcing near 1:3 resonance, for 
different patterns with N modes evenly spaced in Fourier space. Since the energy depends smoothly 
on the angles 0j n (cf. (13)) little change in the energy is expected if the modes are not quite evenly 
spaced. 

For K = 2 and a — — 1 the energies of the patterns with more than 3 modes are very close to each 



other, as seen in Fig. 10 1. For clarity we show therefore in Fig. 10 3 the difference Tn — As one 
increases p, first stripe patterns (N = 1), then square patterns (N = 2), then hexagons (N = 3), and 
eventually supersquares (N = 4) have the lowest energy. For the respective parameter values, these 
are therefore the patterns we expect eventually to arise from noisy initial conditions. Also indicated 



in Figs 10 i,b is the linear stability of these patterns. Here we determine the linear stability of an 
iV-mode pattern by calculating the linear stability of that pattern within the space spanned by the 
N evenly spaced modes and an additional mode rotated an arbitrary angle i(j with respect to the 
ki-mode in Fourier space. For larger p patterns comprised of 5 and 6 modes are linearly stable but 
do not have the lowest energy, though for p > 2 the energies of the patterns with 4 and 5 modes are 
very close. 
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Figure 10: (a) Rescaled energies for evenly spaced modes in the case K = 2. Parameters as in 
FigjlJ) with a = — 1 and $ = 37r/4. Solid (dashed) lines denote linearly stable (unstable) patterns, 
(b) Energy difference Fn — f° r the same data as shown in (a). The 4-mode pattern is preferred 
for p > 1.2. 



In Sec |4.1| we showed that for a more positive a relatively small values of the 1:3- forcing strength 
p are sufficient to stabilize rectangle patterns over a wide range of (cf. Fig. [6^i). As a result we 
expect that the resonant triad can stabilize patterns comprised of more modes than was possible in 
the case a = — 1 depicted above. FigfTTfr shows that this is indeed the case. Increasing a to a — 2.5 
reduces for p > 1 the energy of the 5-mode patterns below that of the 4-mode pattern and renders 
it the pattern with the lowest energy among the patterns with equally spaced modes. We therefore 
expect that in numerical simulations 5-mode patterns would arise over a large range in p. 

For the parameter set that is relevant for the Belousov-Zhabotinsky reaction [48] we found that 
the mode competition is also sufficiently reduced to suggest the stability of multi-mode patterns (cf. 
Figj7|). This is confirmed in Fig 12 , which shows the energies for a = —0.4 and (3 = 1.4. While 
for a = 4 the 4-mode pattern has an energy that is only barely below that of the hexagons and 
reaches those lower values only for large forcing strengths p (Fig 12), increasing a to a = 6 pushes 
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Figure 11: Rescaled energies !Fn for evenly spaced modes in the case K = 2. Parameters as in 
FigjljD with a = 2.5 and & = 37r/4. 5-mode pattern preferred for p > 1. 



the energy of the 4- mode pattern well below that of the hexagons (Figjl2p). We therefore expect 
that the complex patterns we discuss in this paper are accessible in this experimental system. 
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Figure 12: Pattern selection for Belousov-Zhabotinsky parameters [48]. Rescaled energy difference 
J-n — J~a for N evenly spaced modes. Here a = —0.4, (3 = 1.4 , a = 4 (a) and <r = 6 (b). 



For K = 2 cos(tan _1 (2/5)) with $ = 37r/4 (Fig. 13 el), stripes are again stable for small p. As p 
is increased stripes lose stability and square patterns (N = 2) become linearly stable and have the 
lowest energy beyond p w 1.6. Hexagon patterns (TV = 3) become linearly stable at p = 2.7, but 
do not have lower energy than the square patterns. Patterns with more than 3 modes are linearly 
unstable, which is consistent with the predictions based on the coupling coefficient shown in FigjH^i: 
Fourier modes subtending angles smaller than tt/3 compete strongly and therefore do not coexist 
stably. In simulations starting from random initial conditions, we therefore expect either stripe 
patterns (for p < 1.6) or square patterns to arise. For $ = 7r/4, shown in Fig. 13 d, stripes again 
are stable and have the lowest energy for small p. For p = 0.28 they become unstable to rectangle 
patterns spanned by modes subtending an angle of = r = 2 tan -1 (2/5). As p is increased to 
p = 0.46 this rectangle pattern becomes subcritical (cf. Fig{9| and the weakly nonlinear analysis 
taken to cubic order and the associated energy arguments are not sufficient to make predictions 
about pattern selection. For p > 0.28 none of the patterns with equally spaced modes has lower 
energy than the rectangle pattern selected by the resonant triad. Note that in the cubic truncation 
( 13 ) the energy for a pattern diverges when the respective bifurcation changes direction. 



5 Numerical Simulations in Large Domains 

To confirm our predictions for the pattern selection, we perform numerical simulations of the complex 
Ginzburg-Landau equation Being interested in the formation of complex patterns comprised of 
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Figure 13: Rescaled energies Tn for evenly spaced modes in the case K = 2 cos(tan _1 (2/5)). 
Parameters as in FigjTJi with a = -1 and (a) $ = 3tt/4 (cf. Figj8|, (b) $ = tt/4 (cf. Fig{9|. In (a), 
solid (dashed) lines denote linearly stable (unstable) patterns. 



3 or more modes, we focus on the case K — 2. The linear parameters are as in Fig{TJi with nonlinear 
parameters a — —1, $ = 37r/4, and various values of p. We use periodic boundary conditions and 
employ a pseudo-spectral method with Crank-Nicolson-Adams-Bashforth time stepping 

First, to test the weakly nonlinear analysis we focus on the regime where both hexagons and 4- 
mode patterns are linearly stable (FigflQ|) and run numerical simulations in domains of minimal size 
for each pattern. Thus, all participating modes lie exactly on the critical circle. Since the 4-mode 
pattern with evenly spaced modes does not lie on a regular Fourier grid we approximate it by modes 
spaced at X = 2 tan -1 (1/2) w 53° and 2 = tt/2 - 2 tan -1 (1/2) w 37° apart. Due to the smooth 
dependence of the energy on the angles 0j n (cf. Eq.(13)) little change in the energy is expected if 



the modes are not quite evenly spaced. From FigfTofr, for a l:3-forcing strength p = 1 hexagons and 
4-mode patterns are both linearly stable to stripes and the 4-mode patterns are also linearly stable 
to square patterns. The numerical simulations confirm these stability predictions. Starting with a 
noisy initial condition generated from a uniform distribution with amplitude 0.05 we get as expected 
hexagons (Figfl4^i) and supersquares (FigjlZfr). Fig 14 also illustrates the slight diffference in the 
domain sizes used in the simulations, which are required to accomodate the patterns. 



o 



o 




(a) 



b) 



Figure 14: Small-system numerical simulations of (|6| for 1:3 forcing strength p = 1. (a) Hexagon 
pattern in rectangular domain, L x = 47r/k SH « 23.49, L y = L x /\^3 « 13.56. (b) Supersquare 
pattern in square domain, L = 27r/(cos(tan _1 (1/2))^^^^ « 13.13. Parameters as in Figjl]3 (case 
K = 2) with a = -1, $ = 3tt/4. 
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To investigate the competition between iV-mode patterns with different values of N in the same 
computational domain we perform simulations in a large system of linear size L x = L y = L, given 
below, representing 40 wavelengths. FigfTo]o shows that near p = 1 hexagons and 4-mode patterns are 
both linearly stable, but the pattern with minimal energy depends on p. In order to investigate the 
competition between these planforms we start each simulation with the same noisy initial condition, 
generated from a uniform distribution with amplitude 0.01, and vary only p. Moreover, in order to 
clearly identify the nonlinear competition between hexagons and 4-mode patterns and its dependence 
on p, we choose L such that the Fourier modes with = tt/3 and 6 = 7r/4 have the same growth 
rates. This is achieved by taking L = 40(27r /(k^ SH>} — d)) « 473.39 for which the modes for the 
hexagons and the 4-mode patterns are at an equal distance d = 0.004159 on opposite sides of the 
critical circle. 



x 10 4 ^ x 10 4 




Figure 15: Strobed time dependence of the power spectrum in an annulus around the critical circle 
for varying 1:3 forcing strengths p. Parameters as in FigfTj} with a = — 1 and $ = 37r/4. All 
simulations use the same noisy initial condition, (a) p = 0.8, (b) p = 1, (c) p = 1.2, (d) p = 2, (e) 
p = 3. 



To visualize the competition between the modes near the critical circle, Fig |15] shows the temporal 
evolution of the magnitude of the Fourier modes inside a narrow annulus around the critical circle, 
divided into bins with an angular width of 2° plotted stroboscopically at multiples of the period 
Att/v. For p = 0.8, based on energy arguments and the linear stability illustrated in Fig 10 3, we 
anticipate the final state to consist of a stripe pattern. This is confirmed by the numerical simulation. 
While in the Fourier transform Fig 15 1 initially three modes come up, reflecting the linear stability of 
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hexagons for p = 0.8, ultimately the pattern with lower energy, stripes, dominates and only a single 
mode remains. The competition can also be seen in in TimeEvolutionMovie_rhoIs0p8.mov, which 
shows the temporal evolution of the full pattern strobed with period Air/v. More careful inspection 
shows that the middle peak, which seems weaker, actually consists of two modes of lesser strength 
that are slowly converging to a single strong peak. The splitting of the peak into two modes reflects 
a slight undulation of the resulting pattern, which is apparent in the snapshot of part of the final 



solution shown in Fig 16 1. 

Figfl5]3 shows the evolution of the Fourier transform for p = 1 corresponding to the time evolution 
of the full pattern shown in TimeEvolutionMovie_rhoIsl.mov. As with the p = 0.8 case, initially 
the amplitudes of three modes grow. In contrast to the p = 0.8 case, however, for p = 1 stripes are 
unstable to hexagons, and hexagons have the lowest energy. Correspondingly the three modes persist 
resulting in the hexagon patterns shown in Fig (IB}}. Note that the hexagon patterns are comprised 
of domains of up- and down-hexagons characterized by white and black centers, respectively. This 



reflects the fact that the amplitude equations (11 ) for these subharmonic patterns have no quadratic 



term, so neither the up- nor the down-hexagons are preferred. Whether eventually one of the two 
types wins out over the other depends on the interaction between the fronts connecting the domains. 
For p = 1.2 we anticipate from FigfTofc that the resulting solution will be a 4-mode pattern. 



Indeed, Fig|15p shows 4 modes in the Fourier transform, though they are not quite equally spaced. 
The corresponding pattern evolution is shown in TimeEvolutionMovie_rhoIslp2.mov; Fig(T6fc shows 
the final state, characterized by supersquare (dash-dotted, blue circle) and antisquare (dashed, yellow 
circles) elements with approximate 4-fold rotational symmetry [52], as well as approximate 8-fold 
symmetric elements (solid, white circle). 

From FigfTofr a 4-mode pattern is also anticipated for p = 2, but the numerical simulation 
actually results in a 5-mode pattern as shown in Figjl5]i. The energies for evenly-spaced 4- and 
5-mode patterns are very close for this forcing strength. Therefore changes in the energies that 
result from an uneven distribution of the modes may render the 5-mode pattern energetically lower. 
The corresponding pattern features elements with approximately 5- and 10- fold rotational symmetry 



(dashed, yellow and solid, white circles, respectively, in Fig 161). For p = 3 even more modes persist 
for a long time and at the final time shown, Airt/v = 5 • Kr, the pattern has not settled yet into 
its final state. The temporal evolution of the full patterns for p = 2 and p = 3 are shown in 
TimeEvolutionMovie_rhoIs2.mov and TimeEvolutionMovie_rhoIs3.mov, respectively. 




Figure 16: Zoom view (0.25L x 0.25L) of final state at time Airt/u = 5 • 10 4 for varying values of 1:3 
forcing p with L = 4Q(2n / (k^ SH ^ — d)) ~ 473.39. Other parameters as in Figjlj} (case K = 2) with 
a = — 1, <l> = 37r/4. The circles mark characteristic features of the patterns, (a) p = 0.8, (b) p = 1, 
(c) p = 1.2, (d) p = 2, (e) p = 3. 
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To characterize the ordering process taking place in the evolution of the patterns we calculate 
the spectral pattern entropy S = J2 i -pij lnp^-, where p^ = P~ij/(J2ijPij) 18 the normalized power 
in the Fourier mode with amplitude pij and the sum includes all modes within the critical annulus 
binned into 10° segments. The entropy allows an estimate of the effective number of Fourier modes 
e s contributing to the pattern. Fig(l7] illustrates the temporal evolution of e s for different values 
of the 1:3 forcing strength p. The number of significant modes increases with p, although not 
monotonically for all times: for example, the curves corresponding to p = 1 and p = 1.2 cross 
several times, near dimensionless times vt/Air = 12 x 10 3 , 17 x 10 3 ,35 x 10 3 , before beginning to 
converge smoothly to e s ~ 3 and 4, respectively. While the decrease in the effective number of 
Fourier modes with time in the transients is monotonic it occurs in spurts. For p = 3, for instance, 
there is a sudden dip near dimensionless time vt/Ait = 45 x 10 3 . These dips are related to the sudden 



disappearance of modes that is apparent in the time evolution plots Fig 15i-e. 




10 15 20 25 30 35 40 45 50 55 60 
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Figure 17: Effective number of Fourier modes e s plotted as a function of time for different values of 
the 1:3 forcing strength p. Data taken from numerical simulations of (|6| for parameters as in FigfTj} 
(case K = 2) with a = -1, $ = 37r/4. 

The complex Ginzburg-Landau equation (|6| describes the evolution of the oscillation amplitude 
on a slow time. The full time evolution of the underlying system involves also the Hopf frequency. 
To illustrate the full time dependence the movie SimulationOver2Periods_withFastOscs.mov shows 
the evolution of the pattern over two forcing periods obtained for p — 2 as one might see it in an 
experiment. Fig [18] shows the corresponding temporal evolution of the real and imaginary parts 
of one of the spatial Fourier modes. The beating reflects the quasi-periodic forcing of the system. 
Details about the quantity shown in the movie SimulationOver2Periods_withFastOscs.mov and in 



Fig 18 are given in the Appendix. 



6 Conclusion 

We have demonstrated analytically and confirmed numerically that in systems that are near a 
Hopf bifurcation to spatially uniform oscillations forcing with judiciously chosen waveforms can 
stabilize complex periodic and quasi-periodic patterns comprised of up to 4 and 5 Fourier modes. 
Essential for the success of this approach was the use of a quasi-periodic forcing function in which 
the frequency content near twice the Hopf frequency consists of two frequencies; thus, the forcing 
is slowly modulated in time. Since the investigated complex patterns are subharmonic in time with 
respect to that modulation they arise in a pitch-fork bifurcation and are amenable to a weakly 
nonlinear analysis. The stabilization of the complex patterns was achieved by exploiting resonant 
triad interactions that result from the additional forcing component with a frequency close to three 
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Figure 18: Temporal evolution of the real and imaginary parts of one of the spatial Fourier modes 



of the pattern shown in Fig 16 3 and in the movie SimulationOver2Periods_withFastOscs.mov. 



times the Hopf frequency. 

As expected from investigations in the context of the Faraday system [19,24-26,28,53], our 
weakly nonlinear analysis showed that the resonant triad interaction can significantly modify the 
competition between Fourier modes of different orientation through the excitation of weakly damped 
modes. This damping is controlled by the differences in frequency and by the amplitude of the two 
forcing components near twice the Hopf frequency. It can be chosen to enhance the self-coupling 
of the modes significantly, which effectively reduces the competition between modes over a quite 
wide range of angles subtended by those modes. As in the Faraday system, this stabilizes complex 
patterns in the forced oscillatory systems investigated here. Alternatively, the forcing can be chosen 
to reduce the competition between modes differing in their orientation be a quite specific angle r , 
which can be tuned over quite some range. In the Faraday system this approach can also lead to 
complex 4-mode patterns, which effectively consist of the combination of two square patterns that 
are rotated by an angle close to r . In the forced oscillatory systems investigated here this mechanism 
alone does not yield complex patterns since in the unforced system squares are never stable, while 
in the Faraday system square patterns are stable for low viscosity [4,23,24]. Thus, in the present 
context this mechanism only leads to the stabilization of rectangle patterns. For sufficiently strong 
forcing their bifurcation can even be made subcritical, even though the stripe patterns (and the 
Hopf bifurcation itself) are supercritical. 

We have complemented the weakly nonlinear analysis by direct simulations of the extended 
complex Ginzburg-Landau equation in large domains to study the competition between different, 
linearly stable, complex patterns. As expected from the variational character of the amplitude 
equations and the dependence of the energy of the various patterns on the forcing strength, we find 
that as the resonant triad interaction is increased more complex patterns dominate over simpler 
patterns. For the parameters chosen in the numerical simulations we find patterns with 4-fold 
symmetric elements reminiscent of super-squares and anti-squares [52] as well as 5-fold symmetric 
elements. The weakly nonlinear analysis (see Figj5j)) suggests that for other system parameters 
patterns comprised of more modes yet could be stable. 

It should be noted that we reached the regime in which complex patterns are stable by tuning 
the amplitudes and phases of the forcing function, which constitute external forcing parameters. 
Therefore we expect that the complex patterns should be accessible quite generally in forced oscilla- 
tory systems, in particular also in chemical oscillators [37,40,54]. Specifically, our weakly nonlinear 
analysis indicates that the dependence of the patterns on the nonlinear dispersion coefficient a of 
the unforced system can be compensated by the strength of the forcing close to three times the Hopf 
frequency (cf. Figj5}}). Moreover, for the parameters determined experimentally in the oscillatory 
Belousov-Zhabotinsky reaction [48] we have shown explicitly that moderate forcing strengths are suf- 
ficient to stabilize 4-mode patterns. It should be noted, however, that the complex patterns possibly 
arise only very close to onset. Therefore the application of our results to experiments may require 
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systems with relatively large aspect ratios and a very careful tuning of the forcing parameters. 

To characterize the temporal evolution of the patterns starting from random initial conditions 
and to distinguish the resulting patterns quantitatively we used the spectral entropy of the patterns, 
which quantifies the effective number of Fourier modes of the patterns. It allowed a clear distinction 
between patterns with three, four, or more significant modes. However a more detailed quantitative 
characterization of the transients that captures also the competition between multi-mode structures 
like super-squares and anti-squares, which have the same number of participating modes, is still an 
open problem (cf. [55]). The long-time scaling of the ordering process of such complex structures 
and a comparison with the ordering in stripe [56,57] or hexagon patterns [58] may also be interesting 
to study. Most likely, the number of different types of defects and their mutual interaction may play 
an important role. However, for the subharmonic patterns discussed here the amplitude equations 
do not contain any terms of even order. Therefore the strong interaction between defects that is 
associated with those terms and that should make their dynamics in particular interesting in the 
case of 5-mode patterns [59] will not be present here. 

Here our focus was on systems below the Hopf bifurcation (/i < 0). It would be interesting to 
pursue a similar study above the Hopf bifurcation. There the spontaneous oscillations, which do 
not arise in the Faraday system, and their competition with phase-locked patterns driven by the 
forcing may lead to additional complexity. For single-frequency forcing above the Hopf bifurcation, 
labyrinthine stripe patterns arise from the oscillations through front instabilities and stripe nucle- 
ation [54]. It is unknown what happens if the stripes are unstable to the more complex patterns 
discussed here. 

It would also be interesting to consider opposite signs for the detunings U21 and ^22 of the 1 : 2- 
forcing, frequency so that one of the detunings is just above the Hopf frequency, and the other just 
below. Then only one of the two forcing terms will yield a pattern- forming instability [36], while 
the other term will induce a phase-locking of spatially homogeneous oscillations. The result of the 
interaction between these two instabilities is not known. 

Finally, we did not explore the subcritical cases bo < or b(0)/bo < — 1. The weakly nonlinear 
analysis to cubic order is insufficient to make predictions about pattern selection in these cases. 
Fig J9] suggests that for K = 2 cos(tan _1 (2/5)) large-amplitude rectangle patterns with 6 close to 
r = tan -1 (2/5) could be stabilized by the forcing 

We gratefully acknowledge discussions with A. Rucklidge and M. Silber. This work was supported 
by NSF grants DMS-322807 and DMS-0309657. 



Appendix 



In Section [2] we derived the form of the extended CGLE Q based on symmetry arguments. To 
illustrate that the additional terms arise from a general forcing function in a natural way we derive 
the coefficients of the equation here for the Brusselator, which often has served as a simple model 
for chemical oscillations [39]. With the forcing included it can be written as 



du 
~dt 

dv 
dt 



1-(1 + B)u + D u Au + (l + / 2 (cos(x) sin((2 + v 2 i)t) + sin(x) sin((2 + i/ 22 )*))+ (14) 
+/i sin((l + u t )t) + f 3 sin((3 + u 3 )t)) u 2 v, 



Bu + D v Av — u v, 



(15) 



where u and v represent two reacting and diffusing chemical species. In the formulation ( 14|15 ), 
which differs slightly from the form presented in [39] in the coefficients of the linear terms, the Hopf 
bifurcation occurs at B = 2 and the Hopf freque ncy is given by a; = 1. The small parameters 
i>j = 0(5 2 ), S <C 1, represent the detuning (cf. eqs.(3|4)), the parameters j\ = 5 s fi, / 2 = 5/ 2 , and 
f 3 = Sfs represent the small 1:1-, 1:2-, and 1:3 forcing strengths with fj = (9(1), and the parameters 



D u and D v are the diffusion coefficients. To derive the extended CGLE ([6| we expand ( 14|15 ) near 
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the Hopf bifurcation about the basic state (u,v) = (1, B) as 

") = ( b) +S (( <_I T i)/2 )c(i.Me"- + c.c.+ ( <-'-f/" )/ s e3"- + c . c .) +0 (^). 
To obtain the equation for the complex amplitude A we extract the frequency ^21 and write 

C = yf2/3 e i ( tan ~ 1 ( 1 / 8 )+ u ^ t )/ 2 A. (16) 



After rescaling the spatial coordinates as (x,y) = \J (D u + D v )/2(x, y) we arrive at Eq.(|6| with the 
coefficients given by 

fi = V2 + 15|/ 3 | 2 /8, (17) 

a = -(^ 21 /2 + 33|/ 3 | 2 /8), (18) 

7 = 6/ 2 /V65, (19) 

p = I/3IV205/24, (20) 

= tan-^M/^ + arg^). (21) 

As before ^ = ^22 — ^21 and 77 = pe 1 ^ . 

To give an impression of the temporal evolution of the patterns as they may be seen in ex- 
periments we show in (SimulationOver2Periods_withFastOscs.mov) (see also Fig ,18) a movie of the 
^-component of the Brusselator. More precisely, we show only the spatial dependence associated with 
the wavevectors on the critical circle. In terms of the expansion ( p~Q| ) the corresponding amplitudes 
can be written as 



Uk =z (7-1-1)^1 g (X„ + iy„)e-V^/2W 

+ i~\ + vf £ {Xn ~ iY ny nut e- lU2it/2 ~ ibji ) ■ 

^ ' n= — 00 / 

Here Z represents the steady-state amplitude of the Fourier mode given by the fixed-point solu- 
tion of Eq{ll] 
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